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1.  Locations  of  18,101  earthquakes  (red)  and  1770  explosions  (blue)  in  southern 

California  from  2000  to  2005  as  recorded  by  broadband  stations  (yellow)  of  the  Southern 
California  Seismic  Network  (SCSN).  2 

2.  A  cartoon  showing  how  measured  spectra  can  be  modeled  as  a  product  of  event, 

station,  and  travel-time  dependent  terms.  3 

3.  Stacked  P-wave  source  displacement  spectra  from  2000  to  2005  within  bins  of 

estimated  seismic  moment  for  17810  earthquakes  and  1744  quarry  blasts.  (A)  Stacked 
earthquake  source  terms  obtained  from  the  iterative  inversion.  Red  line  shows  the 
empirical  Green’s  function  (EGF)  used  to  correct  these  spectra  for  attenuation  and  other 
path  effects  assuming  a  constant  stress  drop  model.  (B)  EGF  corrected  earthquake  source 
terms  compared  to  predictions  of  the  Madariaga  (1976)  source  model  (dashed  lines).  (C) 
Stacked  source  terms  for  quarries.  4 

4.  (1)  Two  examples  of  earthquake  waveforms  and  spectra.  (2)  Three  examples  of  quarry 

blast  waveforms  and  spectra,  (a)  Waveforms  windowed  around  the  P-wave  first  arrival 
recorded  on  the  vertical  component.  Event  ID  and  station  acronym  are  shown  in  the 
upper  left,  (b)  Waveforms  windowed  around  the  S-wave  arrival  recorded  on  the  rotated 
transverse  component  of  same  station,  (c)  Spectra  for  P  (solid  red),  S  (solid  blue),  and 
respective  noise  levels  (dashed  red  and  dashed  blue),  (d)  EGF-corrected  P-wave  source 
spectra  (red)  together  with  the  best-fitting  source  model  (dashed).  5 

5.  Left:  RMS  Misfit  vs.  corner  frequency  for  earthquakes  (red)  and  quarry  blasts  (blue). 

Right:  RMS  misfit  vs.  seismic  moment.  6 

6.  Moment  vs.  corner  frequency  of  earthquakes  (red)  and  quarry  blasts  (blue).  The  dashed 

lines  represent  stress  drop  estimates  according  to  the  Madariaga  (1976)  model.  7 

7.  Histograms  comparing  the  distribution  of  RMS  misfit  to  an  co‘  source  model  for  both 

earthquakes  (top)  and  quarry  blasts  (bottom).  The  dashed  vertical  lines  divide  the 
distributions  into  10%  and  90%  parts.  Note  that  90%  of  the  quarries  have  model  misfits 
greater  than  that  of  90%  of  the  earthquakes.  8 

8.  Histograms  comparing  the  distribution  of  corner  frequency  for  earthquakes  (top), 

shallow  earthquakes  (middle)  and  quarry  blasts  (bottom).  The  dashed  vertical  lines 
divide  the  distributions  into  10%  and  90%  parts.  9 

9.  Histograms  comparing  the  distributions  of  S/P  amplitude  ratios  between  earthquakes 

(top)  and  quarry  blasts  (bottom).  The  left  panels  include  only  traces  with  an  epicentral 
distance  of  at  least  100  km.  10 
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10.  Locations  of  ten  specific  event  clusters  labeled  as  quarry  blasts  (red).  Each  cluster 
consists  of  a  different  number  of  individual  events.  1 1 


11.  Example  quarry  1  (see  Fig.  10  for  location).  Upper  left:  Stacked  quarry  blast  source 

spectra  in  0.2  local  magnitude  bins.  The  red  line  shows  the  EGF  computed  from 
earthquake  spectra.  Upper  middle:  EGF  corrected  source  spectra.  Upper  right: 
Histograms  of  RMS  misfit  and  corner  frequency  for  all  earthquakes  (red)  and  the  quarry 
example  (blue).  The  bottom  panels  from  left  to  right  show  moment  vs.  corner  frequency, 
RMS  misfit  vs.  corner  frequency,  and  RMS  misfit  vs.  moment  for  all  earthquakes  (red) 
and  the  example  quarry  (blue).  12 

12.  Example  quarry  2  (see  Fig.  10  for  location).  Upper  left:  Stacked  quarry  blast  source 

spectra  in  0.2  local  magnitude  bins.  The  red  line  shows  the  EGF  computed  from 
earthquake  spectra.  Upper  middle:  EGF  corrected  source  spectra.  Upper  right: 
Histograms  of  RMS  misfit  and  corner  frequency  for  all  earthquakes  (red)  and  the  quarry 
example  (blue).  The  bottom  panels  from  left  to  right  show  moment  vs.  corner  frequency, 
RMS  misfit  vs.  corner  frequency,  and  RMS  misfit  vs.  moment  for  all  earthquakes  (red) 
and  the  example  quarry  (blue).  13 

13.  Example  quarry  3  (see  Fig.  10  for  location).  Upper  left:  Stacked  quarry  blast  source 

spectra  in  0.2  local  magnitude  bins.  The  red  line  shows  the  EGF  computed  from 
earthquake  spectra.  Upper  middle:  EGF  corrected  source  spectra.  Upper  right: 
Histograms  of  RMS  misfit  and  corner  frequency  for  all  earthquakes  (red)  and  the  quarry 
example  (blue).  The  bottom  panels  from  left  to  right  show  moment  vs.  corner  frequency, 
RMS  misfit  vs.  corner  frequency,  and  RMS  misfit  vs.  moment  for  all  earthquakes  (red) 
and  the  example  quarry  (blue).  14 

14.  Example  quarry  4  (see  Fig.  10  for  location).  Upper  left:  Stacked  quarry  blast  source 

spectra  in  0.2  local  magnitude  bins.  The  red  line  shows  the  EGF  computed  from 
earthquake  spectra.  Upper  middle:  EGF  corrected  source  spectra.  Upper  right: 
Histograms  of  RMS  misfit  and  corner  frequency  for  all  earthquakes  (red)  and  the  quarry 
example  (blue).  The  bottom  panels  from  left  to  right  show  moment  vs.  corner  frequency, 
RMS  misfit  vs.  corner  frequency,  and  RMS  misfit  vs.  moment  for  all  earthquakes  (red) 
and  the  example  quarry  (blue).  15 

15.  Example  quarry  5  (see  Fig.  10  for  location).  Upper  left:  Stacked  quarry  blast  source 

spectra  in  0.2  local  magnitude  bins.  The  red  line  shows  the  EGF  computed  from 
earthquake  spectra.  Upper  middle:  EGF  corrected  source  spectra.  Upper  right: 
Histograms  of  RMS  misfit  and  corner  frequency  for  all  earthquakes  (red)  and  the  quarry 
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16.  Example  quarry  6  (see  Fig.  10  for  location).  Upper  left:  Stacked  quarry  blast  source 
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Histograms  of  RMS  misfit  and  corner  frequency  for  all  earthquakes  (red)  and  the  quarry 
example  (blue).  The  bottom  panels  from  left  to  right  show  moment  vs.  corner  frequency, 
RMS  misfit  vs.  corner  frequency,  and  RMS  misfit  vs.  moment  for  all  earthquakes  (red) 
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1.  SUMMARY 


Routine  seismic  discrimination  between  earthquakes  and  explosions  has  been  a  long¬ 
standing  goal  in  nuclear  test  ban  treaty  research  (for  a  recent  review,  see  Stump  et  ah, 
2002).  A  variety  of  methods  have  been  employed,  including  amplitude  ratios  among 
regional  phases  (e.g.,  Bennett  and  Murphy,  1986;  Wuster,  1993;  Plafcan  et  ah,  1997; 
McLaughlin  et  ah,  2004),  spectral  studies  (e.g.,  Taylor  et  ah,  1988;  Gitterman  and  van 
Eck,  1993;  Kim  et  al.,  1994;  Walter  et  ah,  1995;  Gittennan  et  al.,  1998),  coda  studies 
(e.g.,  Su  et  al.,  1991;  Hartse  et  al.,  1995),  ripple-fire  detection  schemes  (e.g.,  Hedlin  et  al., 
1990;  Smith,  1993;  Carr  and  Garbin,  1998;  Hedlin,  1998;  Arrowsmith  et  al.,  2006),  and 
other  methods  (e.g.,  Musil  and  Plesinger,  1996;  Parolai  et  al.,  2002;  Leidig  et  al.,  2004; 
Tibuleac  et  al.,  2004). 

The  goal  of  this  project  is  to  systematically  analyze  and  compare  source  spectra  from 
locally  recorded  earthquakes  and  explosions  in  southern  California  (Figure  1)  in  order  to 
develop  new  insights  into  discrimination  methods.  Advances  in  data  storage  and 
computer  capabilities  make  possible  much  more  extensive  analyses  than  have  been 
performed  in  the  past,  which  will  provide  a  better  picture  of  the  distribution  of  source 
spectral  properties  and  amplitudes.  By  examining  tens  of  thousands  of  events,  we  will 
quantitatively  characterize  differences  between  earthquakes  and  explosions  in  terms  of 
their  spectral  content  and  their  P/S  energy  ratios.  We  also  plan  to  identify  and  examine 
anomalous  events,  in  particular  earthquakes  that  may  appear  like  explosions  in  spectral 
discrimination  methods  in  order  to  determine  how  common  they  are  and  whether 
alternate  discrimination  techniques  can  be  applied. 

The  project  builds  upon  a  recently  completed  large-scale  analysis  of  southern 
California  earthquake  spectra  (Shearer  et  al.,  2006),  to  include  a  set  of  1770  mining  and 
other  explosions  between  2000  and  2005.  The  Shearer  et  al.  earthquake  study  has  already 
provided  the  largest  set  of  earthquake  spectra  and  stress  drops  computed  to  date,  showing 
that  individual  event  stress  drops  range  between  0.2  and  20  MPa.  The  large  number  of 
stations  and  events  available  in  southern  California  make  possible  empirical  calibration 
methods  to  remove  receiver  response  and  path  propagation  effects.  Our  efforts  focus  on 
southern  California  because  of  the  unmatched  size  and  quality  of  the  available  data,  but 
we  expect  the  results  and  insights  will  be  applicable  to  other  regions  of  more  direct 
interest  to  nuclear  monitoring  programs.  While  the  Shearer  et  al.  (2006)  study  analyzed 
1989-2001  data  from  short-period  vertical-component  stations,  we  examine  2000-2005 
data  from  three-component,  broadband  stations.  The  newer  data  have  the  advantage  of 
the  horizontal  components  and  a  larger  dynamic  range  (i.e.,  the  older  data  clip  on 
earthquakes  above  -M3.5). 
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Figure  1.  Locations  of  18,101  earthquakes  (red)  and  1770  explosions  (blue)  in 
southern  California  from  2000  to  2005  as  recorded  by  broadband  stations  (yellow)  of  the 
Southern  California  Seismic  Network  (SCSN). 


2.  TECHNICAL  APPROACH 

The  Southern  California  Seismic  Network  (SCSN)  has  several  hundred  stations  and 
records  about  12,000  to  35,000  earthquakes  each  year.  Recently  we  began  storing 
seismograms  from  all  archived  events  in  an  online  RAID  system  that  provides  rapid  and 
random  access  to  the  data  (Hauksson  and  Shearer,  2005).  Spectra  are  computed  as 
follows:  For  each  seismogram  we  pick  the  P  and  S  arrivals  and  estimate  their  amplitudes. 
This  is  done  using  the  operator  pick,  if  available,  or  using  the  output  of  an  automatic 
picking  algorithm  for  a  window  around  the  predicted  arrival  time  (based  on  the  catalog 
event  location  and  a  1-D  velocity  model).  Traces  are  resampled  to  a  uniform  100  Hz 
sample  rate.  Spectra  are  computed  for  1.28  s  noise  and  signal  windows,  immediately 
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before  and  after  the  pick  time.  We  compute  results  for  all  available  channels  and 
components  for  both  P  and  S,  including  rotation  of  the  horizontals  (if  present)  into 
transverse  and  radial  records.  Both  signal  and  pre-event  noise  spectra  are  corrected  to 
displacement  and  stored  in  a  special  binary  format. 


Event  term:  et 


i  Station  term:  s, 


✓'Travel  time  term: 


Figure  2.  A  cartoon  showing  how  measured  spectra  can  be  modeled  as  a  product  of 
event,  station,  and  travel-time  dependent  terms. 


We  apply  a  signal-to-noise  (STN)  cutoff  to  the  spectra,  requiring  that  the  STN 
amplitude  ratio  be  at  least  5  for  three  separate  bands  of  5  to  10  Hz,  10  to  15  Hz  and  15  to 
20  Hz.  Next,  we  process  the  spectra  in  order  to  isolate  source,  receiver  and  propagation 
path  effects.  This  is  an  important  step  because  individual  spectra  tend  to  be  irregular  in 
shape  and  difficult  to  fit  robustly  with  theoretical  models.  However,  by  stacking  and 
analyzing  thousands  of  spectra  it  is  possible  to  obtain  more  consistent  results.  The  basic 
approach  is  illustrated  in  Figure  2  and  is  similar  to  that  used  by  Warren  and  Shearer 
(2000,  2002)  and  Prieto  et  al.  (2004).  Each  observed  displacement  spectrum  dy(f)  from 
source  i  and  receiver  j  is  a  product  of  a  source  term  e,  (which  includes  the  source 
spectrum  and  near-source  attenuation),  a  near-receiver  term  sj  (which  includes  any 
uncorrected  part  of  the  instrument  response,  the  site  response  and  the  near-receiver 
attenuation),  and  a  travel-time  dependent  term  tkaj)  (which  includes  the  effects  of 
geometrical  spreading  and  attenuation  along  the  ray  path).  In  the  log  domain,  this 
product  becomes  a  sum: 

dij  Si  T  Sj  "F  tkfij)  "F  }  jj 

where  ry  is  the  residual  for  path  ij.  We  parameterize  t  in  terms  of  the  predicted  P  travel 
time  between  the  source  and  receiver,  using  the  event  locations  and  velocity  model  from 
Lin  et  al.  (2007).  This  accounts  for  both  the  event  depth  and  the  source-receiver  distance 
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The  travel-time  term  tk(ij)  is  discretized  by  its  index  k  at  1  s  increments  in  travel  time. 
Because  each  station  records  multiple  events  and  each  event  is  recorded  by  multiple 
stations,  this  is  an  over-determined  problem.  We  solve  this  equation  using  a  robust, 
iterative,  least-squares  method  in  which  we  sequentially  solve  individually  for  the  terms 
4,  Sj,  and  eh  keeping  the  other  terms  fixed  at  each  stage.  We  suppress  outliers  by 
assigning  Ll-norm  weights  to  misfit  residuals  greater  than  0.2  s  (or  less  than  -0.2  s).  This 
weighting  scheme  is  necessary  to  ensure  robustness  with  respect  to  a  small  number  of 
spectra  with  large  excursions  compared  to  the  bulk  of  the  data.  In  practice  we  found  that 
the  method  converged  rapidly  to  a  stable  solution  after  a  few  iterations. 

Radiation  pattern  differences  are  not  included  and  would  be  difficult  to  include  in  our 
processing  because  they  are  not  generally  available  for  the  smaller  magnitude  events.  By 
using  multiple  stations  for  each  source,  however,  radiation  pattern  effects  will  tend  to 
average  out.  Note  that  this  method  resolves  only  differences  in  the  relative  shapes  of  the 
spectra.  Without  additional  modeling  assumptions,  it  cannot,  for  example,  resolve  how 
much  of  the  spectral  falloff  is  due  to  source  effects  and  how  much  is  due  to  attenuation 
common  to  all  paths.  The  advantage  of  the  method,  however,  is  that  it  identifies  and 
removes  anomalies  that  are  specific  to  certain  sources  or  receivers.  Because  there  may  be 
difficulties  in  obtaining  reliable  and  accurate  instrument  response  functions  for  many  of 
the  stations  in  the  archive,  this  is  an  important  processing  step  that  provides  a  way  to 
correct  for  some  of  these  problems. 


(a)  (b)  (c) 
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Figure  3.  Stacked  P-wave  source  displacement  spectra  from  2000  to  2005  within  bins  of 
estimated  seismic  moment  for  17810  earthquakes  and  1744  quarry  blasts.  (A)  Stacked 
earthquake  source  terms  obtained  from  the  iterative  inversion.  Red  line  shows  the 
empirical  Green’s  function  (EGF)  used  to  correct  these  spectra  for  attenuation  and  other 
path  effects  assuming  a  constant  stress  drop  model.  (B)  EGF  corrected  earthquake  source 
terms  compared  to  predictions  of  the  Madariaga  (1976)  source  model  (dashed  lines).  (C) 
Stacked  source  terms  for  quarries. 
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Figure  4.  (1)  Two  examples  of  earthquake  waveforms  and  spectra.  (2)  Three  examples  of 
quarry  blast  waveforms  and  spectra,  (a)  Waveforms  windowed  around  the  P-wave  first 
arrival  recorded  on  the  vertical  component.  Event  ID  and  station  acronym  are  shown  in 
the  upper  left,  (b)  Waveforms  windowed  around  the  S-wave  arrival  recorded  on  the 
rotated  transverse  component  of  same  station,  (c)  Spectra  for  P  (solid  red),  S  (solid  blue), 
and  respective  noise  levels  (dashed  red  and  dashed  blue),  (d)  EGF-corrected  P-wave 
source  spectra  (red)  together  with  the  best-fitting  source  model  (dashed). 


3.  RESULTS  AND  DISCUSSION 

Our  focus  has  been  on  the  stacked  source  spectra,  et,  which  we  ultimately  use  to 
estimate  the  moment  and  corner  frequency  of  each  event.  At  this  stage,  however,  the 
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source  spectra  only  contain  relative  information  among  the  different  events.  In  order  to 
estimate  absolute  spectra  from  our  source  stacks,  we  use  the  local  magnitude  Ml  to  obtain 
the  scaling  factor  necessary  to  convert  our  relative  moment  estimates  to  absolute  moment 
and  we  use  an  empirical  Green's  function  approach  to  correct  the  spectral  shapes  for 
attenuation  and  other  path  effects  (for  details,  see  Shearer  et  ah,  2006).  To  study  the 
average  shape  of  the  spectra,  we  stack  our  results  within  equally  spaced  bins  in  estimated 
seismic  moment  (obtained  from  the  low-frequency  part  of  the  spectrum).  Figure  3  shows 
these  stacked  spectra  for  both  earthquakes  and  quarry  blasts.  The  dashed  lines  show  the 
best- fitting  predictions  of  the  of  source  model  of  Madariaga  (1976),  assuming  a  constant 
stress  drop. 
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Figure  5.  Left:  RMS  Misfit  vs.  comer  frequency  for  earthquakes  (red)  and  quarry  blasts 
(blue).  Right:  RMS  misfit  vs.  seismic  moment. 

Figure  3  shows  that  averaged  earthquake  spectra  in  southern  California  are  well  fit  by 
a  standard  source  model.  However,  the  averaged  quarry  spectra  appear  anomalous  in  at 
least  two  respects:  (1)  They  exhibit  large  misfit  compared  to  the  source  model 
predictions,  and  (2)  They  have  generally  steeper  falloffs  at  high  frequencies  than  co‘  , 
which  will  lead  to  lower  corner  frequencies  and  stress  drop  estimates.  The  lack  of  high 
frequency  radiation  from  the  quarries  is  somewhat  surprising  and  may  reflect  ripple  firing 
and/or  strong  near-surface  attenuation.  The  effect  is  also  apparent  in  individual  source 
spectra,  as  illustrated  in  Figure  4.  Here,  we  show  example  waveforms  and  spectra  for 
earthquakes  and  events  labeled  as  quarry  blasts,  at  good  and  fair  signal-to-noise  ratios. 
The  spectra  for  the  two  earthquakes  fit  the  theoretical  source  spectra  much  better  than  the 
quarry  blast  spectra.  We  also  notice  that  the  signal-to-noise  ratio  is  fairly  poor  for  most  S- 
wave  spectra.  We  attribute  this  to  contamination  from  P  coda.  We  have  therefore  focused 
our  study  on  the  P-wave  spectra  obtained  from  vertical  components. 

In  any  case,  we  attempt  to  use  these  two  differences  to  discriminate  between  earthquakes 
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and  quarry  blasts  in  southern  California.  We  do  this  by  computing  the  best- fitting  of 
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source  model  to  the  individual  EGF-corrected  source  spectra.  For  each  event,  we  obtain 
an  estimate  of  the  moment,  the  comer  frequency  and  a  measure  of  the  RMS  misfit  to  the 
source  model.  Figure  5  shows  the  spectral  misfit  plotted  against  the  corner  frequency  and 
the  seismic  moment.  First  of  all,  we  notice  that  the  spectral  misfit  does  not  depend  on 
either  the  comer  frequency  of  the  spectmm  or  the  seismic  moment  of  the  event,  which 
makes  this  parameter  viable  to  be  used  as  a  discriminant.  Note  that  the  quarry  blasts  have 
generally  higher  misfits  and  smaller  comer  frequencies  than  the  explosions.  We  also 
observe  that  earthquakes  span  a  much  wider  moment  range  than  the  man-made 
explosions.  In  general,  the  two  populations  are  not  completely  separated  and  there  is 
some  degree  of  overlap,  particularly  in  the  comer  frequency  estimates.  Figure  6  shows 
the  separation  of  the  two  populations  in  the  corner  frequency-moment  domain,  with  lines 
of  constant  stress  drop  indicated.  We  obtain  generally  lower  stress  drop  estimates  for  the 
quarry  blasts,  even  though  stress  drop  is  not  defined  for  an  artificial  explosion.  We  also 
observe  that  the  distribution  of  earthquake  stress  drops  implies  self-similarity  (stress  drop 
not  dependent  on  moment),  although  the  scatter  increases  toward  smaller  moments. 
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Figure  6.  Moment  vs.  corner  frequency  of  earthquakes  (red)  and  quarry  blasts  (blue).  The 
dashed  lines  represent  stress  drop  estimates  according  to  the  Madariaga  (1976)  model. 


We  next  look  at  histograms  of  the  earthquake  and  quarry  blast  distributions  with  respect 
to  RMS  misfit  (Fig.  7)  to  see  how  well  the  populations  separate.  We  find  an  empirical 
value  of  0.2  RMS  misfit  at  which  the  earthquake  and  quarry  blast  populations  separate  at 
the  90%  level.  In  other  words,  90%  of  the  earthquakes  show  an  RMS  misfit  of  <  0.2, 
whereas  90%  of  the  quarry  blasts  show  an  RMS  misfit  of  >  0.2,  A  similar  histogram 
with  respect  to  comer  frequency  (Fig.  8)  shows  an  overall  weaker  separation,  where  the 
10%  and  90%  quantile  cannot  be  drawn  at  the  same  corner  frequency  value.  In  order  to 
mle  out  the  possibility  that  the  lower  corner  frequency  observed  for  quarry  blasts  is  an 
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effect  of  the  increased  attenuation  in  the  shallow  subsurface  and  thus  attributable  to  event 
depth,  we  compare  the  distribution  of  corner  frequencies  with  shallow  earthquakes  only 
(Fig.  8  b).  We  note  that  there  is  no  significant  difference  in  the  corner  frequency 
distribution  if  the  earthquake  population  is  limited  to  depths  shallower  than  5  km.  This 
confirms  that  the  lower  corner  frequencies  observed  for  quarry  blasts  are  a  source-related 
effect  rather  than  an  effect  of  near-source  attenuation. 
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Figure  7.  Histograms  comparing  the  distribution  of  RMS  misfit  to  an  co'“  source  model 
for  both  earthquakes  (top)  and  quarry  blasts  (bottom).  The  dashed  vertical  lines  divide 
the  distributions  into  10%  and  90%  parts.  Note  that  90%  of  the  quarries  have  model 
misfits  greater  than  that  of  90%  of  the  earthquakes. 


Our  results  from  analysis  of  S-wavc  spectra  from  transverse-component  records  have  so 
far  been  inconclusive,  in  part  because  of  the  generally  lower  signal-to-noise  ratios  at  high 
frequencies  for  the  S  waves  compared  to  the  P  waves.  However,  the  amplitude  ratio 
between  P-  and  ^-arrivals  is  an  often-used  discriminant  for  explosions.  The  idea  here  is 
that  explosion  sources  preferentially  excite  compressional  waves,  whereas  earthquakes 
are  mainly  due  to  shear  on  a  fault,  and  thus  radiate  S  waves  in  higher  proportion. 


In  order  to  investigate  whether  the  S/P  ratio  could  be  used  as  a  discriminant  in  southern 
California,  we  pick  the  peak  amplitudes  of  the  P-  and  S-arrivals  on  seismograms  with  a 
signal-to-noise  ratio  of  at  least  three  for  the  P  arrival.  This  analysis  is  conducted  on 
seismograms  that  are  filtered  between  1  and  10  Hz.  We  obtain  the  P  pick  from  the 
vertical  component  and  the  S  pick  from  the  rotated  transverse  component. 
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Figure  8:  Histograms  comparing  the  distribution  of  corner  frequency  for  earthquakes 
(top),  shallow  earthquakes  (middle)  and  quarry  blasts  (bottom).  The  dashed  vertical  lines 
divide  the  distributions  into  10%  and  90%  parts. 

Plotting  the  amplitude  ratios  as  a  histogram,  separately  for  quarry  blasts  and  natural 
earthquakes  (Fig.  9),  we  observe  slightly  lower  average  S/P  amplitude  ratios  for  the 
quarries  using  simple  peak  amplitudes.  However,  the  separation  between  the  two 
populations  is  not  apparent.  We  try  to  quantify  the  difference  between  the  two 
populations  by  fitting  a  Gaussian  distribution  to  the  histogram,  from  which  we  obtain  a 
mean  and  standard  deviation  for  the  P/S  ratios  of  natural  earthquakes  and  man-made 
explosions. 

We  observe  that,  for  southern  California,  quarry  blasts  and  earthquakes  are  not 
distinguishable  by  their  S/P  amplitude  ratios.  One  reason  for  this  could  be  the 
interference  of  P  coda  into  the  S  arrivals.  In  order  to  avoid  this  possibility,  we  look  at 
only  records  with  an  epicentral  distance  of  greater  than  100  km  (see  Fig.  10).  This 
ensures  that  P  and  S  arrivals  are  sufficiently  separated  in  time  in  order  that  the  S  arrival  is 
not  contaminated  by  P  coda.  This  amplitude  ratio  is  comparable  with  a  Pn/Sn  ratio.  We 
notice  that  the  separation  between  the  two  populations  is  much  more  pronounced  when 
we  restrict  the  amplitude  ratio  analysis  to  epicentral  distances  greater  than  100  km, 
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roughly  equivalent  to  the  critical  distance  for  the  Pn  phase.  We  observe  that  at  these 
distances,  the  S/P  ratios  decrease  for  the  quarry  blasts,  whereas  they  increase  for  the 
earthquakes.  It  should  be  noted,  however,  that  restricting  the  epicentral  distance  to  greater 
than  100  km  could  result  in  a  bias  because  we  are  in  effect  restricting  out  database  to  look 
only  at  larger  events.  Overall,  we  notice  that  the  S/P  amplitude  ratios  are  not  very  useful 
as  a  discriminant  between  earthquakes  and  explosions  in  southern  California. 


logioS/P  ratio 


Figure  9.  Histograms  comparing  the  distributions  of  S/P  amplitude  ratios  between 
earthquakes  (top)  and  quarry  blasts  (bottom).  The  left  panels  include  only  traces  with  an 
epicentral  distance  of  at  least  100  km. 


Next,  we  look  at  individual  clusters  of  quarry  blast  events  that  can  be  attributed  to 
individual  operations.  Most,  but  not  all,  of  these  locations  can  be  attributed  to  mining.  We 
look  at  10  individual  event  clusters  that  are  highlighted  in  Figure  10.  Cluster  1  through  5, 
as  well  as  cluster  No.  8  are  located  in  the  Mojave  Desert  and  can  be  correlated  with 
surface  mining  operations  visible  in  satellite  imagery.  Clusters  6,7  and  10  are  located  in 
densely  populated  areas  of  Riverside  and  San  Diego  County  and  cannot  be  easily 
correlated  with  a  particular  operation  or  surface  feature.  Some  of  these  events  may  be  due 
to  grading  and/or  construction.  Cluster  No.  9  is  more  distributed  in  nature  and  covers  the 
area  of  the  China  Lake  Naval  Air  Weapons  Station.  We  suspect  that  testing  of  military 
equipment  could  be  adding  to  the  catalog  in  this  region.  However,  this  is  also  an  area 
with  a  high  rate  of  natural  seismicity,  so  it  is  likely  that  some  events  in  this  region  may  be 
mislabeled.  Figures  1 1  to  20  represent  summary  plots  for  each  individual  cluster, 
showing  the  stacked  spectra  per  magnitude  bin,  their  respective  fit  to  the  standard  source 
model,  as  well  as  histograms  and  cross  plots  for  the  misfit  and  the  comer  frequency.  As 
before,  we  observe  the  generally  steeper  falloff  rate  of  the  spectra,  leading  to  a  high 
misfit  with  respect  to  an  earthquake  source  model.  These  plots  also  help  to  identify 
possible  characteristic  patterns  in  a  particular  cluster.  Such  patterns  could  be  indicative  of 
particular  blasting  practices  (i.e.,  ripple-fired  explosions  with  characteristic  delay  times, 
surface  vs.  underground  explosions,  etc.).  Such  characteristic  patterns  are  not  directly 
evident  from  the  plots.  However,  we  do  note  slight  differences  between  some  of  the 
clusters.  For  example,  cluster  8  shows,  on  average,  larger  moment  releases  than,  e.g., 
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cluster  3.  This  could  be  an  indication  either  of  different  charge  sizes  used  in  the  two 
mining  operations,  or  it  could  be  related  to  the  ground  coupling  (i.e.,  hard  vs.  soft  rock). 
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Figure  10.  Locations  of  ten  specific  event  clusters  labeled  as  quarry  blasts  (red).  Each 
cluster  consists  of  a  different  number  of  individual  events. 
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Figure  11.  Example  quarry  1  (see  Fig.  10  for  location).  Upper  left:  Stacked  quarry  blast 
source  spectra  in  0.2  local  magnitude  bins.  The  red  line  shows  the  EGF  computed  from 
earthquake  spectra.  Upper  middle:  EGF  corrected  source  spectra.  Upper  right: 
Histograms  of  rms-misfit  and  corner  frequency  for  all  earthquakes  (red)  and  the  quarry 
example  (blue).  The  bottom  panels  from  left  to  right  show  moment  vs.  corner  frequency, 
rms-misfit  vs.  corner  frequency,  and  rms-misfit  vs.  moment  for  all  earthquakes  (red)  and 
the  example  quarry  (blue). 
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Figure  12.  Example  quarry  2  (see  Fig.  10  for  location).  Upper  left:  Stacked  quarry  blast 
source  spectra  in  0.2  local  magnitude  bins.  The  red  line  shows  the  EGF  computed  from 
earthquake  spectra.  Upper  middle:  EGF  corrected  source  spectra.  Upper  right: 
Histograms  of  rms-misfit  and  corner  frequency  for  all  earthquakes  (red)  and  the  quarry 
example  (blue).  The  bottom  panels  from  left  to  right  show  moment  vs.  corner  frequency, 
rms-misfit  vs.  corner  frequency,  and  rms-misfit  vs.  moment  for  all  earthquakes  (red)  and 
the  example  quarry  (blue). 
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Figure  13.  Example  quarry  3  (see  Fig.  10  for  location).  Upper  left:  Stacked  quarry  blast 
source  spectra  in  0.2  local  magnitude  bins.  The  red  line  shows  the  EGF  computed  from 
earthquake  spectra.  Upper  middle:  EGF  corrected  source  spectra.  Upper  right: 
Histograms  of  rms-misfit  and  corner  frequency  for  all  earthquakes  (red)  and  the  quarry 
example  (blue).  The  bottom  panels  from  left  to  right  show  moment  vs.  corner  frequency, 
rms-misfit  vs.  corner  frequency,  and  rms-misfit  vs.  moment  for  all  earthquakes  (red)  and 
the  example  quarry  (blue). 
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Figure  14.  Example  quarry  4  (see  Fig.  10  for  location).  Upper  left:  Stacked  quarry  blast 
source  spectra  in  0.2  local  magnitude  bins.  The  red  line  shows  the  EGF  computed  from 
earthquake  spectra.  Upper  middle:  EGF  corrected  source  spectra.  Upper  right: 
Histograms  of  rms-misfit  and  corner  frequency  for  all  earthquakes  (red)  and  the  quarry 
example  (blue).  The  bottom  panels  from  left  to  right  show  moment  vs.  corner  frequency, 
rms-misfit  vs.  corner  frequency,  and  rms-misfit  vs.  moment  for  all  earthquakes  (red)  and 
the  example  quarry  (blue). 
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Figure  15.  Example  quarry  5  (see  Fig.  10  for  location).  Upper  left:  Stacked  quarry  blast 
source  spectra  in  0.2  local  magnitude  bins.  The  red  line  shows  the  EGF  computed  from 
earthquake  spectra.  Upper  middle:  EGF  corrected  source  spectra.  Upper  right: 
Histograms  of  rms-misfit  and  corner  frequency  for  all  earthquakes  (red)  and  the  quarry 
example  (blue).  The  bottom  panels  from  left  to  right  show  moment  vs.  corner  frequency, 
rms-misfit  vs.  corner  frequency,  and  rms-misfit  vs.  moment  for  all  earthquakes  (red)  and 
the  example  quarry  (blue). 
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Figure  16.  Example  quarry  6  (see  Fig.  10  for  location).  Upper  left:  Stacked  quarry  blast 
source  spectra  in  0.2  local  magnitude  bins.  The  red  line  shows  the  EGF  computed  from 
earthquake  spectra.  Upper  middle:  EGF  corrected  source  spectra.  Upper  right: 
Histograms  of  rms-misfit  and  corner  frequency  for  all  earthquakes  (red)  and  the  quarry 
example  (blue).  The  bottom  panels  from  left  to  right  show  moment  vs.  corner  frequency, 
rms-misfit  vs.  corner  frequency,  and  rms-misfit  vs.  moment  for  all  earthquakes  (red)  and 
the  example  quarry  (blue). 
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Figure  17.  Example  quarry  7  (see  Fig.  10  for  location).  Upper  left:  Stacked  quarry  blast 
source  spectra  in  0.2  local  magnitude  bins.  The  red  line  shows  the  EGF  computed  from 
earthquake  spectra.  Upper  middle:  EGF  corrected  source  spectra.  Upper  right: 
Histograms  of  rms-misfit  and  corner  frequency  for  all  earthquakes  (red)  and  the  quarry 
example  (blue).  The  bottom  panels  from  left  to  right  show  moment  vs.  corner  frequency, 
rms-misfit  vs.  corner  frequency,  and  rms-misfit  vs.  moment  for  all  earthquakes  (red)  and 
the  example  quarry  (blue). 
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Figure  18.  Example  quarry  8  (see  Fig.  10  for  location).  Upper  left:  Stacked  quarry  blast 
source  spectra  in  0.2  local  magnitude  bins.  The  red  line  shows  the  EGF  computed  from 
earthquake  spectra.  Upper  middle:  EGF  corrected  source  spectra.  Upper  right: 
Histograms  of  rms-misfit  and  corner  frequency  for  all  earthquakes  (red)  and  the  quarry 
example  (blue).  The  bottom  panels  from  left  to  right  show  moment  vs.  corner  frequency, 
rms-misfit  vs.  corner  frequency,  and  rms-misfit  vs.  moment  for  all  earthquakes  (red)  and 
the  example  quarry  (blue). 
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Figure  19.  Example  quarry  9  (see  Fig.  10  for  location).  Upper  left:  Stacked  quarry  blast 
source  spectra  in  0.2  local  magnitude  bins.  The  red  line  shows  the  EGF  computed  from 
earthquake  spectra.  Upper  middle:  EGF  corrected  source  spectra.  Upper  right: 
Histograms  of  rms-misfit  and  corner  frequency  for  all  earthquakes  (red)  and  the  quarry 
example  (blue).  The  bottom  panels  from  left  to  right  show  moment  vs.  corner  frequency, 
rms-misfit  vs.  corner  frequency,  and  rms-misfit  vs.  moment  for  all  earthquakes  (red)  and 
the  example  quarry  (blue). 
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Quarryl  0 


Figure  20.  Example  quarry  9  (see  Fig.  10  for  location).  Upper  left:  Stacked  quarry  blast 
source  spectra  in  0.2  local  magnitude  bins.  The  red  line  shows  the  EGF  computed  from 
earthquake  spectra.  Upper  middle:  EGF  corrected  source  spectra.  Upper  right: 
Histograms  of  rms-misfit  and  corner  frequency  for  all  earthquakes  (red)  and  the  quarry 
example  (blue).  The  bottom  panels  from  left  to  right  show  moment  vs.  corner  frequency, 
rms-misfit  vs.  corner  frequency,  and  rms-misfit  vs.  moment  for  all  earthquakes  (red)  and 
the  example  quarry  (blue). 
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4.  CONCLUSIONS 


Earthquake  and  explosions  in  southern  California  exhibit  significant  differences  in 
their  average  E-wave  spectral  properties.  Quarry  blast  spectra  are  not  well-fit  by  standard 
source  models  and  typically  have  lower  corner  frequencies  and  anomalously  steep 
falloffs  at  high  frequencies  compared  to  earthquakes  of  the  same  estimated  moment.  We 
can  therefore  establish  the  RMS  misfit  between  theoretical  spectra  calculated  for  a  Brune- 
type  source  model  and  the  actually  observed  source  spectra  as  a  discriminant  between 
earthquakes  and  explosions  in  Southern  California.  However,  we  were  unable  to 
discriminate  unambiguously  between  earthquakes  and  explosions.  In  particular,  the  two 
populations  still  overlap  to  such  an  extent  that  not  only  an  earthquake  could  be 
misclassified  as  an  explosions,  but  also  (which  is  potentially  worse)  and  explosion  could 
remain  undetected  by  being  misclassified  as  an  earthquake.  On  the  other  hand,  we  have  to 
assume  that  the  dataset  used  in  this  analysis  (SCSN  catalog)  did  already  include  a  number 
of  misclassified  events,  since  the  flagging  of  quarry  blasts  in  the  SCSN  database  is  based 
only  on  event  location  and  the  daytime/nighttime  distribution.  It  is  in  this  respect  not  the 
ideal  dataset  to  investigate  an  unambiguous  discrimination  unless  a  more  rigorous 
flagging  of  events  can  be  achieved.  Future  results  from  analysis  of  .S'-wavc  spectra  may 
provide  additional  discriminants. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


EGF 

RMS 

SCSN 

STN 


Empirical  Green’s  Function 
Root  Mean  Square 

Southern  California  Seismic  Network 
Signal-To-Noise 
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